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Field and computer techniques for stem analysis adaptable to both young- and 
old-growth conifers are presented. Field instructions include a step-by-step 
explanation of tree cutting, sectioning, ring count, and measurement techniques. A 
computer program adapted to field techniques reproportions a tree's radial measure- 
ment data, calculates height-age-site index information, punches the reproportioned 
and height-age-site index information on cards, and plots height-age and stem profile 


graphs. 


800 years old and to either American or metric measures. 


These stem analysis techniques are adaptable to trees of any size up to 


These instructions 


enhance the usefulness of the stem analysis research method in obtaining growth 


information for forest managers. 


KEYWORDS: Stem analysis, age determination, ring measurement, data recording 


methods, old-growth conifers. 
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Stem analysis studies within the upper- 
slope mixed-conifer forests of the Cascade 
Range in Oregon and Washington were begun 
in 1965. Objective of these studies was to 
develop site quality and comparative growth 
information among the several coniferous 
species growing in the true fir-hemlock 
forests. Work was done within the natural 
range of noble fir between Stevens Pass 
in Washington and McKenzie Pass in 
Oregon. Stem analysis was selected as a 
technique for developing productivity in- 
formation because it permitted use of 
overmature trees to obtain information 
about growth of both young and old trees. 


During the course of the studies, 
efficient methods of data collection and 
use were devised and applied. From that 
experience, detailed instructions of stem 
analysis field methods, including tree 
felling, sectioning, and measurement of 
sequential radial growth (SRG), are pre- 
sented. Three different methods of 
collecting multiple-ring SRG data are 


explained. One method is used as an 
example to give the reader a complete set 
of instructions for field collection of data 
compatible for automatic data processing 
(ADP). Modification of both data collection 
and computer programing would be neces- 
sary for the other two methods. 


Many old and new stem analysis 
techniques were applied in appropriate 
succession to develop efficient operations. 
The goal of these stem analysis techniques 
is to code and record data compatible for 
ADP. Output for that ADP consists of 
reproportioned SRG measurements, stem 
analysis graph (stem profile chart), and 
height-age graphs and listings for individual 
trees. Any or all of the output items may 
be selected on a single run. 


Both field and computer techniques are 
adaptable to either the American or metric 
system of measurement--only minor 
change in programing is necessary for 
conversion. 
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Computer-produced stem profile chart for noble fir 
with a breast-height age of 432 years. 
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Introduction 


To develop site quality and other growth 
information essential for proper management 
of the upper-slope true fir-hemlock forests 
of the Cascade Range in Oregon and Wash- 
ington, the U.S. Forest Service and Oregon 
State University began stem analysis 
studies in 1965. This work was aimed 
especially at collection of growth data 
from old-growth noble firl/ for estimation 
of site index (DeMars et al. 1970, Herman 
and DeMars 1970). Noble fir stem 
analysis growth data also were used for 
comparison of methods of site index 
estimation (Curtis et al. 1974a); and 
associated Douglas-fir data were used 
in development of site index information, 
specifically for upper-slope Douglas-fir 
in the Cascade Range (Curtis et al. 1974b). 


Species composition within the upper- 
slope mixed-conifer forests of the Cascade 
Range is complex. Radiolongitudinal 
data produced from these stem analysis 
studies are providing the first information 
about measured growth relationships 
among surviving old-growth dominant 
and codominant individual trees within 
and between defined habitats. With such 
a wide representation of different species-- 
as many as 9 or 10 on each acre--consid- 
erable latitude is available for selection 
of a particular combination for managed 
stands. Comparative growth information 
has proved valuable in other mixed 
forest types (Deitschman and Green 1965, 
Green and Alley 1967, Carmean and 
Vasilevsky 1971) and will prove equally 
valuable in the Pacific Northwest. 


Besides providing comparative growth 
information specifically for mixed coni- 
fers in the Cascade Range of Oregon and 
Washington, field and computer techniques 


1/ See table 2, appendix I, page 22, for 
scientific names of trees. 


presented in this paper are applicable to 
stem analysis of both young- and old-growth 
conifers regardless of where they grow. 
Furthermore, such techniques are adapt- 
able to either the American or metric 
system of measurement--only minor 
change in programing would be necessary 
for conversion. 


Natural forest stands containing noble 
fir are found in the Pacific Northwest 
Cascade Range roughly between Stevens 
Pass in northern Washington and McKenzie 
Pass in central Oregon (fig. 1). The 
area south of McKenzie Pass in the 
Willamette National Forest was excluded 
from the study to avoid confusing growth 
variation caused by hybridization with 
Shasta red fir. However, a few Shasta 
red fir and associated tree species were 
sectioned for comparative studies in 
southwestern Oregon. Table 2, appendix I, 
presents a list of the important upper-slope 
coniferous tree species with their assigned 
numerical computer codes. 


Because few upper-slope forest stands 
of 150 years of age and younger were 
found in the Cascades, older trees were 
selected for study. Our stem analysis 
technique permitted us to use overmature 
trees to obtain information on growth of 
both young and old trees. Furthermore, 
stem analysis as a technique was adapted 


to this growth-evaluation research because 
no uniform and random range of age 


classes existed throughout the noble fir 
range. Except for infrequent, scattered 
young-growth stands arising from some 
catastrophic event such as fire or 
windstorm, most stands were 200 to 400 
years old. 


The stem analysis approach has the 
inherent disadvantage of dependence 
among successive measurements on the 
same sample tree and presupposes that 
selected trees were always dominant 


Figure 1.--Stem analysis study 
locations within the natural 
range of noble fir in Oregon 
and Washington. 


throughout their lives (Schlich 1895, 
Dahms 1963). This disadvantage may 
or may not be important. During our 
quest for upper-slope noble fir and 
Douglas-fir site index information, 
occasional checks of heights of two to 
four dominant noble fir trees at selected 
ages at given locations showed no signif- 
icant changes in dominance. 


There are several advantages to the 
use of the stem analysis method. These 
have been discussed at length by several 
workers over many years. Curtis (1964) 
refers to much of that literature and 
points favorably toward the use of the 
method. 


The disadvantage of dependence among 
successive measurements on the same 
tree can be minor compared with the 
disadvantages associated with the empir- 
ical method of attempting to span a site 
quality interval for a given species by 
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selecting for measurement a great many 
representative temporary forest plots. 
Stem analysis provides a solution to the 
problem of unbalanced site selection over 
a wide range of ages because a uniform 
number of height observations always is 
available for any given age for all selected 
sites. 


Even with such widespread species as 
Douglas-fir (McArdle et al. 1961) and 
western hemlock (Meyer 1937, Barnes 
1962), research workers using stand 
data were not able to select equal numbers 
of stands for all ages over an equal number 
of potential growing sites. Beck and 
Trousdell (1973) discussed and demon- 
strated the effect of that problem of 
disproportional sampling of site and age. 


Procedures for stem analysis of trees 
have been both briefly and comprehensively 
described in many basic and advanced 
forestry measurement textbooks and 


bulletins at least as early as 1895 (Schlich). 
A complete description of stem analysis 
methods and uses in early American 
forestry literature was presented by 
Mlodziansky (1898). Detailed field and 
office procedures included use of such 
data to determine past and current 
volumes of individual trees. Many other 
references to early American forestry 
stem analysis methods and uses are given 
by Spurr (1952) and Turnbull. 2/ Later 
references to the technique and one of 

the most complete contemporary descrip- 
tions of procedures are contained in a 
mensuration text by Husch (1963). 


Our paper describes for the first time 
detailed and complete procedures, includ- 
ing collection of field data and use of that 
data, for stem profiles and height-age 
curves produced by Automatic Data 
Processing (ADP). 


Objectives of this paper are: 


1. To describe the field and laboratory 
procedures actually used in stem 
analyses. 


2. To recommend field and laboratory 
procedures applicable to similar 
studies and forestry education exercises. 


3. To present a stem analysis computer 
program with instructions for its use. 


Techniques of field, laboratory, and 
computer operations presented in this 
paper are purposely detailed to furnish 
step-by-step instructions for those who 
wish to use our stem analysis procedures 
to develop site index information and 
secure stem profile data for growth com- 
parisons and volume determinations. 


2/Kenneth J. Turnbull. Stem analysis tech- 
niques and applications and some Studies of 
second-growth Douglas-fir in western Washington. 
Unpublished M.F. thesis, University of Washing- 
ton, Seattle, 1958. 


Stem Analysis Field Methods 


Hundreds of true fir-hemlock stands 
within the natural range of noble fir were 
visited. From these stands, study locations 
were Selected. Habitat type and understory 
vegetation differentials?. 4/ at varying 
elevations and latitudes were bases for 
selection. Frequently, site continuums 
were sSelected--as from ridgetop to valley 
bottom. 


TREE SELECTION AND PREFELLING 
MARKING AND MEASUREMENTS 


After each location was determined, 
the tallest dominant or codominant tree 
of each species (table 2, appendix I) 
present within approximately one-fourth 
acre was selected for stem analysis. 
Trees with obvious defects such as thin 
crowns, stag, forked, or broken tops, 
and bent and crooked stems were avoided, 
as were severely leaning or unbalanced 
trees. Rot indicators such as conks, 
catfaces, frost cracks, and abnormal stem 
swelling caused trees to be rejected. Any 
tree suspected of decay was "bored" at 
stump height with the rounded tip of the 
chain saw. Expelled wood fragments were 
monitored for rot. Trees with "brooms" 
in the crown, whether caused by dwarf 
mistletoe or some other disease, were 
taken if nothing else was available. 


After trees were selected and before 
they were felled, they were circumscribed 
at stump and breast heights with spray 


3/ Jerry F. Franklin. Study plan for develop- 
ment of a habitat type classification of the true 
fir-hemlock forests of the southern Washington 
Cascades. 51p. Unpublished report, Forestry 
Sciences Laboratory, Corvallis, Oregon, 1962. 


ei Jerry Forest Franklin. Vegetation and 
soils in the subalpine forests of the southern 
Washington Cascade Range. Ph. D. thesis, 
Washington State University, Pullman, 132 p. 
1966. 


paint, in planes at right angles to the 
vertical axes of stems. Stump and breast 
height diameters were measured with a 
diameter tape. Where information about 
stem cross-section eccentricity in re- 
lation to cardinal direction was desired, 
special stump and breast height stem 
marks were made before trees were 
felled. North and south (or east and 
west) cardinal direction stem exposures 
were marked with vertical spray paint 
lines (fig. 2). Analyses of cross-section 
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eccentricity information are made 
separately from the same ring count and 
ring interval measurement data collected 
for height-diameter and height-age growth 
evaluations. 


Cardinal direction paint mark--at 180° 
opposite this mark is another 
vertical paint line marked “S” 


TREE FELLING 


Where possible, felling was done along 
contour--not over drainages or rock 
outcrops--to avoid excessive breakage. 
Felling into adjacent trees was avoided as 
much as possible to prevent hangups, severe 
breakage, and loss of candidate treetop. 
Stands with hemlocks were especially 
difficult to work within because the long 
crowns easily became entangled. 


Figure 2.--Diagram show- 
ing painted stump ring, 
breast high ring, and 
cardinal direction, or 
exposure reference 
marks. 
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The stump was cut either 2 inches above 
or 2 inches below the stump ring mark, 
with care taken not to undercut across the 
stump mark. If the cut was below stump 
mark, the undercut left a flat surface on 


base of upper stem section. If cut was 
above stump mark, the undercut left a 
flat surface on stump. The stump section 
was prepared for measurement by cutting 
through the painted ring. When slightly 
leaning or unbalanced trees necessarily 
were selected, special care in stump 
preparation was taken to avoid pulling 
long splinters of wood from the stem or 
stump. Back cut was always maintained 
parallel to undercut. Side cuts were 
made to prevent longitudinal "barber 
chair" splitting of the stem. 


TREE STEM REASSEMBLY 
AND MEASUREMENT 


The tree was limbed from stump to 
top, including all broken stem parts. 
Before length measurements were taken, 
all broken stem pieces were matched to 
provide accurate information on stem 
length. Broken bark and wood patterns, 
wood grain patterns, and knot arrange- 
ments provided ''keys" to accurate 
"reassembly.'' Broken pieces quite often 
were impossible to physically join to- 
gether because of weight and displacement 
distances. Pressurized spray paint cans 
with a fine nozzle provided a good method 
of marking matching points on broken 
stems. 


TREE SECTIONING 


All points on tree stem where section- 
ing cuts were to be made were marked 
with paint. Cuts were made at short log 
length intervals so they could be salvaged. 
Although actual measurements varied, 
the most satisfactory sectioning sequence 
for both salvage and good stem profile 
data was: stump height, 4.5 feet (1.37 
meters) above ground (breast height); 

9 feet (2.74 meters) above b.h.; then 

18 feet (5.47 meters); 36 feet (10. 94 
meters); etc., by 18-foot (5.47-meter) 
increments to a 12-inch (30. 5-centimeter) 
top. Sectioning above the 12-inch top 


was done at shorter invervals--6 to 10 
feet (1.8 to 3 meters)--preferably 
centered near selected internodes. Where 
salvage of stem was not required, section 
cuts above breast height were made at 
individual discretion, but at no greater 
intervals than where salvage was planned. 
In addition to the sectioning points, cuts 
made at a point just below the base of the 
primary live crown in an area undistorted 
by branch swell provide information on 
stem form below and within the crown. 


Section cuts were marked by use of a 
100-foot steel tape graduated in feet and 
hundredths of feet. Average stump height 
was measured and zero end of the tape 
was placed at base of cut stem (allowing 
for the width of the stump section if 
taken from base of stem). Cumulative 
measurements of stem length were 
carried as shown in figure 3. Several 
methods of measuring length were 
tried. All were discarded in favor of 
the system shown in figure 3. Best 
accuracy and, of course, fewer rechecks 
of down trees were necessary when the 
cumulative-height-above-stump measure- 
ments were obtained before the stem 
was bucked and sectioned. Bucking of 
the stem béfore stem length and diam- 
eter measurements are taken is best 
avoided. 
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Figure 3.--Sample of field notes taken in waterproof field surveying book for 
tree measured before it was sectioned. 


Explanation of Field Notes (fig. 3)5/ 


1. Plot name--from some local geographic 


feature or manmade structure such as 
Lower Cedar River, Humpback Moun- 
tain, Yellowjacket Creek, Carpenter 
Lookout Road, etc. 

2. District and Forest. 

3. Plot number. 

4, Species and tree number. 

5. Tree selection date. 

6. Tree felling and sectioning (bucking) 
date. 

7. Stump height above average ground 
level--stump cut was made 1.5 feet 
above ground on uphill side. However, 
recorded stump height above average 
ground level was at least equal to 
1.5 feet. Usually it was greater be- 
cause of sloping ground. 

8. Cumulative stem length from zero 
(stump cut). 

9. Tape ''zeroed" at sectioning point, 
93.4 feet above stump. If tape longer 
than 100 feet is used, ''zeroing tape" 
usually is not required unless stem 
is broken, 


5/ 


= No section at 9 feet above breast height was 
taken on this example because high quality of stem 
dictated salvage for peeler logs. 


Figure 4.--Six bark thick- 
ness measurements were 
taken at approximately 
60° intervals around 
stem circumference on 
each section cut or 
break where d.o.b. 
could not be obtained 
because of bark damage. 
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Cumulative stem length from tape 

zero at 93.4-foot section. 
"Break"'--stem is broken at 32.0" 

from ''93.4" tape zero point. Tape 

is ''zeroed" again and cumulative 
measures continued to next stem 

"break, "' etc. 

Diameter outside bark (d.o. b.) for 
stump and breast height obtained 

before tree was felled, and upper 

stem diameters obtained where bark 

was undamaged before felled stem was 
bucked. 

Where bark was damaged at planned 
point of sectioning, d.o.b., where 
circled, was calculated; double of 
average of six bark thickness measure- 
ments was added to the peeled diameter 
of section. Measurements of bark thick- 
ness were taken at approximately 60° 
intervals on the adjacent, undamaged 
section cut and/or break. Each bark 
thickness measurement represents the 
thickest bark within that 60° arc on the 
stem at or near the section point. Where 
end of cut stem was used, a 6-inch rule 
was used to project undamaged bark sur- 
face to cut stem surface for measurement 
(fig. 4). 


Ruler used to project undamaged 
bark thickness to cut surface of section 
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Diameters inside bark (d.i.b.), 
except where circled, were obtained 
by "slipping" bark from stem section. 
This is easily done during the early 
part of the growing season and soon 
after felling. Later in the growing 
season, all bark thickness was 
measured as described in point 13 
because bark could not be easily and 
completely removed. 

Diameters inside bark, where circled, 
were calculated as described in point 
13 from double of average of six bark 
thickness measurements subtracted 
from d.o.b. 

Total stem length above stump. 

Plot elevation from barometric 
altimeter. 

Average slope of plot area in percent. 
Aspect expressed in quadrant or 
half-quadrant. 

Clear length above stump--no limbs 
or knot indicators on bark to point 
specified. Small epicormic branches 
1/4-inch or less in diameter were 
ignored in determination of clear 
length. 

Length to base of live crown above 
stump--this is to base of primary or 
main live crown exclusive of all green 
epicormic or dead branches. 

Length to base of live crown above 
stump as defined by beginning of 
epicormic branches often found on 
Engelmann spruce and occasionally 
Douglas-fir. 

Crown position--''Dominant, " 
"Codominant," etc. 

Total tree height above average 
ground level. Each tree selected and 
felled for stem analysis provided its 
own individual measurement problems. 
Not all of the previously described 
measurement techniques necessarily 
were used on each tree. In practice, 
the most efficient combination of 
measurement techniques that provided 
accurate data was selected and applied 
as needed. 


DETAILED SECTIONING 
MEASUREMENTS, RING COUNTS, 
AND SEQUENTIAL RADIAL 
GROWTH MEASUREMENTS 


1. Before actual sectioning:9/ 

a. Diameter outside bark was measured 
at all designated sectioning points. 
Where bark was damaged or lost, di- 
ameters outside bark were taken at 
equal distances both ways along log 
from the section cutting mark and 
averaged (except at stump and breast 
height). When this distance was 
greater than twice the diameter at 
that point, the section d.o.b. was 
obtained by d.i.b. measurement and 
individual bark measurements from 
logs adjacent to section. (See explana- 
tory information relative to fig. 4.) 


b. Cardinal direction upper-stem exposure 


marks (see fig. 2). Whichever felled 
stem exposure faces skyward--north, 
east, south, or west--dictates position 
of exposure reference marks across 
all felled stem section points before 
stem is bucked and sectioned. [If tree 
is bucked before exposure marks are 
made, this information quite likely 
will be impossible to collect. 


2. Sectioning: 
Cross sections (disks) 2 to 3 inches 
thick were cut at all marked stem 
section points where diameter was. 
larger than 20 inches (fig. 5). Disks 


6/ If faller worked when measurement crew 
could not be present, many measurement prob- 
lems arose. Trees were felled and bucked to 
length and transverse sections were cut before 
outside bark diameter and cumulative height 
Measurements were taken. Inaccuracies in 
diameter and length measurements taken by 
untrained personnel negated efficiency otherwise 
derived. However, this type of assembly line 
production can be used if felling and bucking 
crews are carefully trained. 


Figure 5.--Cutting stem 
section disk from 
old-growth Douglas-fir. 


1 to 2 inches thick were taken from 
smaller diameter stem sections. Bark 
was removed completely for d.i.b. 
measurement. When only a partial section 
was available because of breakage or buck- 
ing difficulties, its diameter was estimated 
from other associated measurements. A 
partial section was only acceptable when 

it contained the pith and a complete radius 
with an acceptable or ''representative" 
(reasonably normal) ring-growth pattern. 


Trees having annual ring pattern 
destroyed by rot, of course, were rejected. 
However, where decay was encountered 
that had not completely obliterated the 
annual rings, thicker sections of 4 to 6 
inches were cut. Decayed sections were 
carefully handled to prevent breakage. 
Sandwiching them between strips of 
1/4-inch plywood eliminated breakage 
and loss of information. These later were 
carefully soaked in water and frozen to 
facilitate ring counts and measurements 


(Herman et al. 1972a, 1972b). Douglass 
(1928) described an efficient field technique 
which used a paraffin-gasoline solution 

to render decayed wood workable. His 
technique permitted immediate ring 

counts and sequential radial growth 
measurements. 


3. Recommended ring counts and sequen- 
tial radial growth measurements (fig. 6): 


If single ring counts and measurements 
are to provide data for stem profiles, 
data collection is relatively uncompli- 


cated. However, if decadal, half-decadal, 


or other multiple ring measurements 
are to be taken, counts and sequential 
radial growth measurements must be 
taken in a specific manner to provide 
type of data desired. 


At least three different approaches to 
multiple ring counts are possible: 
a. If desirable to compare tree growth 
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Figure 6.--Ring counting using large reading glass. 


of several trees at selected years 

in fixed sequence, multiple-ring count 
must begin at each section circum- 
ference and proceed toward the pith. 
After ring count is completed and ring 
intervals are marked for each section 
radius, sequential radial increment 
measurements are made from the pith 
and recorded on stem analysis ADP 
coding form (fig. 11, appendix J). 

For example, for a tree cut in 1973, 
decadal count and measurement would 
provide growth comparison at all 
section points for 1963, 1953, 1943, 
and so forth back through time at 
10-year intervals. 

b. If desirable to compare tree growth 
at even chronological decades, decadal 
ring count also must begin at section 
circumference and proceed toward 
the pith. However, modification of 
procedure is necessary. For example, 
for trees cut in 1973, if tree growth 


comparison is desirable at even 
decades, i.e., 1970, 1960, 1950, 
etc., the first multiple-ring count 
on each section from circumference 
must span three growth rings to 
reach an even decade.”/ For re- 
mainder of rings to pith, decadal 
counts are taken. Sequential radial 
increment measurements are again 
made from the pith outward along 
each established radius. 

c. If desirable to compare growth of 
trees at specific breast height ages, 
e.g., at 10-year intervals, a more 
complicated procedure of ring count 
and measurement (see footnote 2, p. 3) 
is necessary than has been described. 
Decadal counts and measurements 
for the breast height sections must 
begin at the pzth. For example, 
where decadal ring measurements 


u/ If the 1973 growth ring is present. 


are to be taken, a ring number less than 
10 will likely occur as a partial decadal 
count next to the circumference. The 
number nine or less will be used to begin 
ring counts from the circumference on 

all other sections including the stump. 
Sequential radial increment measurements 
again begin at the pith and continue toward 
the circumference for each radius. 


For this paper, the first-described 
multiple-ring count method--all counts 
begin at the circumference of each 
section--is used. The computer program 


described herein was developed accordingly. 


Some modification of the program is neces- 
sary to produce alternate stem analysis 
information describing growth at even 
chronological decades or at specific 
breast-high age intervals. It should be 
carefully noted that regardless of selected 
multiple ring count approach, all sequential 
growth measurement information is ob- 
tained from cumulative multiple annual 
ring measurements always taken from 

the pith outward toward the circumference 
of each section. For our computer pro- 
gram, that procedure is necessary to 
provide growth data for a normal stem 
profile graph of height over radius. 


Unless the tree was extremely out of 
round or the pith obviously was excessively 
off geometric center, as is sometimes 
found in trees on very steep slopes, only 
one selected radius was taken for sections 
where longest radius differed from 
shortest by 10 percent or less. Otherwise, 
at least two radii--the longest and one 
about 180° to the longest--were selected. 
Three radii, the longest and two others 
about 120° each side of the longest, were 
taken on stump sections. For our stem 
analysis studies, that system always 
resulted in acceptable profiles. Failure 
to select the proper radii for annual ring 
counts and measurements will result in 
unacceptable stem profiles (fig. 15, 
appendix I). 


Successful utilization of stem analysis 
field and laboratory measurements to 
produce stem profiles depends on rea- 
sonably accurate average radii. There 
are several ways to determine the average 
radius (pith to perimeter) of a tree stem 
cross-section or disk. Because of limited 
finances and time, selection of ''repre- 
sentative'' radii was expedient. Most of 
our trees were old-growth specimens, 
and the common occurrence of annual 
ring growth aberrations caused by occluded 
rot, scars, knots, pitch pockets, etc., 
usually prevented actual designation of 
acceptable "true average radii"8/ from 
pith to perimeter on stem section disks. 
Instead, radii were selected that were 
aberration-free and representative 
(ocularly) of the radial growth at each 
section point on the stem. For easy 
measurement of sequential radial growth 
from pith to perimeter and avoidance of 
inaccuracies caused by missing rings 
(Douglass 1928), the shorter radii from 
any section were avoided. 


When more than one radius per section 
was measured, sequential increment 
measurements of each radius were repro- 
portioned (table 1) and a simple arithmetic 
average of each reproportioned radial 
increment was calculated and coded for 


8/ Radii derived from halving diameter tape 
measurements and accepting the small positive 
bias inherent with the diameter tape method 
(Matérn 1956). 
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keypunching (appendix I, fig. 12). Repro- 
portioning was done by the equation: 


Dib/2 


RRIM = TAR 


(RIM) 


Where RRIM = Reproportioned sequential 
radial increment measure- 
ment 

= Diameter inside bark of 
section as determined by 
diameter tape 

= Total length of actual 
measured radius (repre- 
sentative radius) 

= Sequential increment 
measurement from repre- 
sentative radius. 


Dtb 


TAR 


RIM 


Multiple radii were averaged with a 
desk calculator or by a simple program on 
on a desk-top computer. 


A more complicated method of averag- 
ing with the quadratic mean approach was 


suggested by Matérn (1956). According 

to Matérn, the quadratic mean of several 
random radii at each section point pro- 
vides an unbiased estimate of any area, 
convex or not, and is applicable to closures 
such as those projected by tree stem 
cross-sections. We used available repro- 
portioned multiple radius data, and our 
resulting average radii differed little from 
those produced by the quadratic mean 
approach. The differences were small 
(usually second and third decimal place) 
because all measured radii were repro- 
portioned to the section's average radius 
before they were averaged (table 1). 


Reproportionment of sequential radial 
growth measurements for sections with 
a single radius is automatically completed 
by the computer program. Therefore, 
such sections can be coded and keypunched 
as measured. 


Table 1.--Multtple radius computations for WF1, plot 765 


Section: Stump 


Single bark thickness: 1.30 inches 


Section height above ground: 1.5 feet Section ring count: 437 years 
Section? d-o.bi.: 


74.90 inches 


Measured radial data 


Radius 


Radius| Radius| Radius 


1 2 3 
0.78 0.76 0.83 0.73 
1182 1.69 1.80 ae Wal 
Za), 2.74 2.90 2.80 
4.41 3.99 4.19 4.15 
5.96 5.38 5.65 5.61 
36.69 32.65 34.45 34.54 
Si] oS 33.04 34.88 34.97 
37.62 33.46 34.32 35.42 
Sep 33). 72 35.59 35.69 
38.40 34.15 36.05 36.15 
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a PN FR © 


Reproportioned radial 
data 


Radius 


Average 
reproportioned 
Radius Tee whes 
data 
3 

0.83 Oa73 
1.80 Ie 
Zao ots 
4.20 4.19 
Sis 017 5.66 
34.55 34.55 
34.98 34.98 
35.42 35.42 
35.69 35.69 
S6nrS £06 15 


Production of stem profiles from 
nonreproportioned measurement data for 
sections with a single radius, together 
with those having averages of repro- 
portioned multiple radii, is completed 
by the computer program described and 
presented in appendixes II and III. 

Should automatic reproportionment of all 
radii, both single and multiple, be de- 
sired, a revision of that program is 
required. 


RADIAL MEASUREMENTS, 
SEQUENTIAL RADIAL GROWTH 
DATA RECORDS, AND 
SECTION STORAGE 


Where only one radius is required, 
no field counts or measurements need be 
taken. A selected or representative radius 
from a transverse section cut with a small 
chain saw can be taken from the field for 
later processing in the office or labora- 
tory (fig. 7). For sections less than 10 


Figure 7.--Large 
transverse disk 
from sectioned 
noble fir with 
the representative 
radius removed. 


inches in diameter, the complete disk 
should be collected because of difficulty 
in cutting a radial section. 


For both permanent ring count and 
sequential radial growth records where 
complete disk sections cannot be taken, 
counts and measurements on sections 
requiring multiple radii must be completed 
in the field. Heavy clear plastic overlays 
(fig. 8) have been used successfully for 
permanent ring count and growth records. 
Douglass (1919) used a similar technique 
to obtain sequential radial growth records. 
He used "paper rubbings" for transcrib- 
ing ring-growth information from trans- 
verse radial sections. 


Prior to field use, fine parallel lines 
spaced one-half inch apart were drawn 
on the transparent plastic with a sharp 
round point such as a push pin or shortened 
dissection needle. After any radius on a 
section was established, counts were made 
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as previously described and decadal or 
other multiple ring marks were made on 
the wood with a hard sharp-pointed blue 
lead pencil. A moist surface provided 
darker and sharper blue pencil marks. 
For accuracy and uniformity, we recom- 
mend that all marks be made so that the 
edge of the blue pencil line coincides with 
the outside edge of the summer wood for 
each marked ring. Transcribing these 
marks to the fine parallel lines on a 
plastic overlay was relatively simple, 
again done with a sharp-pointed scriber. 
Where any radius was longer than a 
single line on the overlay, a second line, 
etc., with proper documentation was used 
for continuation. These plastic overlays 
were set aside for later checking and 
measurement and are available for other 
observations. Scribed marks on the 
clear plastic are easily seen with a 
dark-colored background under the plastic. 


A transverse section of wood was also 
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“ Figure 8.--Tracing radial 
ring pattern on plastic 
overlay. 


collected, cut with a small chain saw; 

it contained the best representative radius 
from sections (stump, b.h., etc.) from 
which multiple radii were transcribed to 
plastic. These radial transverse sections 
were valuable for correcting any errors 
later encountered. 


Identification of sections is extremely 
important. Aluminum flashing squares 
(such as those used for referencing and 
boundary marking during surveying)-- 
with identification scribed by ballpoint 
pen--were found to be the most durable 
method of marking wood sections. They 
were firmly stapled adjacent to and not 
covering selected radius of each section 
top. Each section was identified with 
plot number, tree species, tree number, 
and stem section location; e.g., 

"P70 NF1 127.3" means plot 70, noble 
fir 1, 127.3 feet above stump. ''Top" of 
each section was always identified. 


Radial transverse sections and disks 
for each tree were packaged with fiber- 
glass tape. An identification tag bearing 
plot, tree number, species, and general 
location was attached to each package 
of radial transverse sections (fig. 9). 


Stem Analysis Laboratory 
and 
Computer Methods 


Annual-ring counts and sequential 
radial growth measurements not com- 
pleted in the field were made in the 
laboratory by the count plan selected. 
The sequential radial growth marks 
transcribed in the field to the plastic 
overlays also were measured according 
to the selected count plan, and the data 
were entered on a coding form (fig. 11, 
appendix I). Data on the coding forms 
were keypunched on standard 80-column 
ADP cards. 


Sequential radial growth was measured 


Figure 9.--Radial trans- 
verse sections and disks 
packaged for laboratory 
examination. Tag identi- 
fies package of sections 
of white fir, tree 
number 1, from the Rogue 
River National Forest. 


with a 24-inch machinist rule graduated 
to the nearest hundredth of an inch except 
where annual rings were closely spaced. 
For transverse sections with rings very 
close together, a dendrochronograph 
equipped with a traveling binocular zoom 
microscope was uSed (fig. 10), and 
sequential radial growth was measured 
to the nearest hundredth of a millimeter. 
To correspond with ruler-measured data 
for our studies, metric information was 
converted to inches. 


For easier ring counts and measure- 
ments along representative radii, the 
wood surface on stem section disks or 
radial transverse sections was smoothed 
with a thin-bladed, chisel-edged, hobby- 
type razor-knife or was sanded with an 
electric sander. Sanding with very fine 
sandpaper usually was done only when 
narrow growth rings (about 40 or more 
per inch) were encountered. Douglass 


(1928) favored the razor-knife method 
of wood surface preparation over that of 
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sanding, and our work supported that 
opinion. To make the annual rings more 
distinct, we used water for most species. 
For resinous wood of pines and occasion- 
ally Douglas-fir, we found that kerosene 
or other light oil was the best agent for 
enhancing the annual rings. Sophisticated 
wood-staining techniques were unnecessary. 


Where there was concern for radial 
shrinkage in the collected radial sections 
and transverse disks, the current year's 
or "green" radii were placed on the plastic 
overlays in the field. Two marks--pith 
center and section perimeter edge--were 
transcribed to these plastic overlays. 
Before laboratory ring counts and sequen- 
tial radial growth measurements were 
taken, green dimensions were restored 
by soaking the collected specimens in 
water (Herman and DeMars 1970) for at 
least 24 hours, then draining them 
overnight. Radii of green sections that 
were correctly marked on the plastic 
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Figure 10.--A dendro- 
chronograph was used 
to measure radial 
growth where annual 
rings were closely 
spaced. 


always corresponded closely to those of 
the water-restored sections. Ring 
intervals for these field-selected repre- 
sentative radii then were marked off 
and measured according to the selected 
count plan. 


COMPUTER PROGRAM FOR 
STEM ANALYSIS DATA 


At least two previous programs using 
stem analysis data for computing and 
graphing tree growth are in existence 
(Brace and Magar 1968, Pluth and 
Cameron 1971). Both of these automated 
methods machine-plot derived tree growth 
parameters. Neither program provides 
a standard height-radius stem profile as 
provided in this paper. Both programs 
provide volume and basal area computa- 
tions that could be adapted for use with 
our program. 


Preparation for computer production 
of stem analysis data requires coding of 
the field and the laboratory data from 
stem disks and scribed plastic overlays. 
Figure 12, appendix I, is a stem analysis 
data code completed for "NF1, plot 75." 
From data sheets such as these, 80-column 
ADP cards were keypunched and verified. 
A computer was programed to produce for 
each tree: (1) a stem analysis graph 
(fig. 13, appendix I) using reproportioned 
radial data and straight-line interpolated 
heights, (2) listing of reproportioned 
radial data corresponding to the stem 
analysis profile, (3) a set of ADP cards 
punched to contain the reproportioned 
data, (4) a height-age listing by the 
selected ring-count interval using 
straight-line interpolation, (5) a height- 
age graph of the original tree data, and 
(6) a set of ADP cards punched to con- 
tain the selected ring count interval 
height-age data. 


Any or all of the previous items can 
be selected as output for an individual 
tree on a Single run. These options 
are explained in appendix II, card 
type I. 


EDITING STEM ANALYSIS DATA 
USING THE COMPUTER PROGRAM 


In actual practice, usually only the 
stem profile was plotted first to dis- 
cover if any errors existed in the raw 
data. Such errors (figs. 14 and 15, 
appendix I) were readily determined 
and corrections were made to insure 
that the changes were correct. Follow- 
ing the production of an acceptable stem 
profile graph, the remaining computer 
program options were selected for 
completion. Data editing using the 
stem profile graph should be done 
cautiously because not all apparently 
eccentric profile graph lines are 
measurement errors. Some line aberra- 
tions can be caused by growth problems 
such as old stem breaks, forks, etc. 


Oe ea ot 


A copy of the stem analysis program 
listing is in appendix II]. Further infor- 
mation about the program may be ob- 
tained from the Pacific Northwest Forest 
and Range Experiment Station, P.O. Box 
3141, Portland, Oregon 97208. 


Literature Cited 


Barnes, George H. 


1962. Yield of even-aged stands of western hemlock. U.S. Dep. Agric. Tech. 


Bull 12/35) oles tllus. 


Beck, Donald E., and Kenneth B. Trousdell 
1973. Site index: accuracy of prediction. 


USDA For. Serv. Res. Pap. SE-108, 


7 p., illus. Southeast. For. Exp. Stn., Asheville, N.C. 


Brace, L. G., and K. M. Magar 


1968. Automated computation and plotting of stem analysis data. Can. Dep. For. & 


Rural Dev. Publ. 1209, 8 p., illus. 


Carmean, Willard H., and Alexander Vasilevsky 


1971. Site index comparisons for tree species in northern Minnesota. USDA For. Serv. 
Res. Pap. NC-65, 8 p., illus. North Cent. For. Exp. Stn., St. Paul, Minn. 


17 


Curtis, Robert O. 
1964. A stem analysis approach to site index curves. For. Sci. 10: 241-256, illus. 


Curtis, Robert O., Donald J. DeMars, and Francis R. Herman 
1974a. Which dependent variable in site index-height-age regressions? For. Sci. 20: 
74-87, illus. 


Curtis, Robert O., Francis R. Herman, and Donald J. DeMars 
1974b. Height growth and site index estimates for Douglas-fir in high-elevation forests 
of the Oregon- Washington Cascades. For. Sci. 20:307-316, illus. 


Dahms, Walter G. 
1963. Correction for a possible bias in developing site index curves from sectioned 
tree data. J. For. 61: 25-27, illus. 


Deitschman, Glenn H., and Alan W. Green 
1965. Relations between western white pine site index and tree height of several 
associated species. USDA For. Serv. Res. Note INT-22, 28 p., illus. Intermt. 
For. & Range Exp. Stn., Ogden, Utah. 


DeMars, Donald J., Francis R. Herman, and John F. Bell 
1970. Preliminary site index curves for noble fir from stem analysis data. USDA 
For. Serv. Res. Note PNW-119, 9p., illus. Pac. Northwest For. & Range 
Exp. Stn., Portland, Oreg. 


Douglass, A. E. 
1919. Climatic cycles and tree growth. A study of the annual rings of trees in relation 
to climate and solar activity. Carnegie Inst. Wash. Publ. No. 289(1), 127 p., illus. 


Douglass, A. E. 
1928. Climatic cycles and tree growth. A study of the annual rings of trees in relation 
to climate and solar activity. Carnegie Inst. Wash. Publ. No. 289(II), 166 p., illus. 


Green, Alan W., and Jack R. Alley 
1967. Ryaluatine species alternatives for National Forest land Sugabie of growing 
western white pine. USDA For. Serv. Res. Pap. INT-43, 40 p., illus. Intermt. 
For. & Range Exp. Stn., Ogden, Utah. 


Herman, Francis R., and Donald J. DeMars 
1970. Techniques and problems of stem analysis of old-growth conifers in the 
Oregon- Washington Cascade Range. Jn J. H. G. Smith and J. Worrall (eds.), 
Tree-ring analysis with special reference to Northwest America, p. 74-77, illus. 
Univ. B. C. Fac. For. Bull. 7. Vancouver. 


Herman, Francis R., Clark E. Smith, and John E. Firth 
1972a. Freezing decayed wood to facilitate ring counts and width measurements. USDA 
For. Serv. Res. Note PNW-187, 4p., illus. Pac. Northwest For. & Range Exp. 
Stn., Portland, Oreg. 


18 


Herman Francis R., Clark E. Smith, and John E. Firth 
1972b. Ring count of decayed wood made easier through freezing. J. For. 70: 743, illus. 


Husch, Bertram 
1963. Forest mensuration and statistics. 474 p., illus. New York: The Ronald Press Co. 


McArdle, Richard E., Walter H. Meyer, and Donald Bruce 
1961. The yield of Douglas-fir in the Pacific Northwest (slightly revised). U.S. Dep. 
Agric. Tech. Bull. 201, 74 p., illus. 


Matérn, Bertil 
1956. On the geometry of the cross-section of a stem. Medd. Fran Statens 
Skogsforskningsinst. Band 46. NR11:1-27, illus. 


Meyer, Walter H. 
1937. Yield of even-aged stands of Sitka spruce and western hemlock. U.S. Dep. Agric. 
Tech. Bull. 544, 86p., illus. 


Mlodziansky, A. K. 
1898. Measuring the forest crop. U.S. Dep. Agric. Div. For. Bull. 20, 71p., illus. 


Pluth, D. J., and D. R. Cameron 
1971. Announcing FORTRAN IV program for computing and graphing tree growth 
parameters from stém analysis. For. Sci. 17: 102. 


Sehlich, William 
1895. A manual of forestry. Vol. III. Forest management. 397 p., illus. London: 


Bradbury, Agnew and Co. 


Spurr, Stephen H. 
1952. Forest inventory. 476 p., illus. New York: The Ronald Press Co. 


19 


a ; 


ban 
ve 
Du 


i 
oy 


i 


th 
yu, 
i 


' paket ® 
ry nnts 
: siete 
7 A 
‘; i . 
1‘ + i 
‘2 
. 
t 
i éV 4 ak 
j 4 
if AF Y : 
a) 1 Sy Re i me 
s ; noe 
Rin) ! Rope Wan iene (hen ok y 


Mi Hie & 


} Has ee at) i 
q ¢ 1 
5 = 
o ’ 
is 


Appendix | 


21 


Table 2.--Noble fir and associated spectes felled and sectioned during stem 
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Common name 


noble fir 2/ 
Douglas-fir— 
Paci fac: sal liver far, 
western hemlock 
mountain hemlock 
western Naee pine 
grand fir 

western redcedar 
subalpine fir 
western larch 
lodgepole pine 
Engelmann spruce 
Shasta red fir 
ponderosa pine2/ 


1/ 


analysts studtes 


: eens 1 
Scientific name 


Species 
ADP code 


Abtes procera Rehd. 022 
Pseudotsuga menztestt (Mirb.) Franco 202 
Abtes amabilis (Dougl.) Forbes 011 
Tsuga heterophylla (Raf.) Sarg. 263 
Tsuga mertenstana (Bong.) Carr. 264 
Pinus monttcola Dougl. 3/119 
Abtes grandis (Dougl.) Lindl. — 015 
Thuja plicata Donn. 254 
Abtes lastocarpa (Hook.) Nutt. 019 
Larix oecetdentalis Nutt. 073 
Pinus contorta Dougl. 108 
Picea engelmannti Parry 093 
Abtes magnifica var. shastensts Lemn. 021 
Pinus ponderosa Laws. 122 


— Scientific names follow Elbert L. Little, Jr., Check list of native 
and naturalized trees of the United States (including Alaska). U.S. Depart- 
ment of Agriculture, Agriculture Handbook 41, 472 p., 1953. 


2/ 


— Species associated with Shasta red fir felled and sectioned south of 
the known normal range of noble fir. Douglas-fir and grand fir (white fir) 
were sampled in both the noble fir range and Shasta red fir range. 


3/ 


— Grand fir and white fir considered as single species for this phase 
of study--both coded ''015." 
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Figure 1l1.--Stem analysis data code sheet. 
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Figure 13.--Computer produced stem profile chart for noble fir 1) DLOE Hor. 
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Figure 14.--Computer produced stem profile graph showing error in data. 
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Figure 15.--Computer produced stem profile showing results of inadequate data 
measurements at lower stem section points. 
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Card Input Formats 


With the exception of the ''00'' card, the program!/ data input cards can be punched 


directly from the stem analysis data code sheets (figs. 11 and 12 of appendix I). 


Descrip- 


tive plot and miscellaneous tree information to be punched on the ''00"' card comes from 


various places in the original field notes. 


following format is mandatory. 


Card Type I--''00" Card 


Implied 
Card decimal 
columns designation 
1-3 XxX 
4 xX 
5-7 XxX 
8-11 XXX. X 
12-15 XXXX 
16-17 XX 
18-19 xXx 
20 xX 
21-22 xx 
23-24 XxX 
25-28 XXX. X 
29-32 XXX. X 
33-36 XXX, X 
37-40 XXX. X 
1/ 


Identification 


Plot number 

Tree number 

Species code 

D.b.h. (o.b. ) 

Breast height ring 
count 

Card number 

Number of sections 
taken for the tree 

Ring count measure- 
ment interval (years) 


Year cut 
Crown class 


Height to first live 
limb (feet) 

Height to live crown 
(feet) 

Clear length (feet) 


Crown length (feet) 


Transcription of that information into the 


Special 


notes 


Table 2, appendix I 
Inches 


Always "00" 
(Including tip, 
limit--30) 
1= 10 yr 
5= Syr 
3= lyr 


Dominant and tallest 
on plot--11 
Dominant--1 
Codominant--2 
Intermediate--3 
Suppressed--4 
Above average 
ground level 
Above average 
ground level 
Above average 
ground level 
Above average 
ground level 


— Regardless of tree size or age, program is adaptable to data from 30 section points 


where the 30th section point is the tip having a diameter of ''0'"' and a single bark thickness 


of ''0.'' Also, program is limited to 80 incremental measurement intervals. 


Expressed 


in terms of stump age, program will accommodate an 80-year-old tree at 1-ring measure- 
ment intervals, a 400-year-old tree at 5-ring measurement intervals, or an 800-year-old 


tree at 10-ring measurement intervals. 
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Card Type I--''00" Card, continued 


Card 
columns 


41-45 
46 
47 


48-50 
51-56 


57-62 
63-65 


66-69 
70-74 
75 


76 


CU 


78 


19 


80 


2/ 


Implied 
decimal 


designation 


XXXXX 
xX 
x 


Identification 


Elevation (feet) 
Blank 
Aspect 


Slope (percent) 
Latitude north 


Longitude west 


Preliminary noble fir 
site index 

Total tree height (feet) 

Blank 

Card output of 
reproportioned data 


Printed output of 
reproportioned data 


Stem profile 


Card output of height-age 
data by the selected ring 
count interval 

Printed output of height-age 
data by the selected ring 
count interval 

Height-age plot of original 
tree data 


— Donald J. DeMars, Francis R. Herman, and John F. Bell. 


curves for noble fir from stem analysis data. 


9 p., illus. 


Special 
notes 


Aspect codes 
Level - 0 
North, northeast 1,2 
East, southeast 3,4 
South, southwest 5,6 
West, northwest 7,8 


Cols. 51-53 degrees 
Cols. 54-56 min 
Cols. 57-59 degrees 
Cols. 60-62 min 
From PNW-119 
Res. Note 19702/ 


To suppress output, 
punch a "1"; otherwise 
leave blank. 

To suppress output, 
punch a "1"; otherwise 
leave blank. 

To suppress output, 
punch a "1"; otherwise 
leave blank. 

To suppress output, 
punch a "1"; otherwise 
leave blank. 

To suppress output, 
punch a "'1"; otherwise 
leave blank. 

To suppress output, 
punch a "1"; otherwise 
leave blank. 


Preliminary site index 


USDA Forest Service Research Note PNW-119, 


Pacific Northwest Forest and Range Experiment Station, Portland, Oregon, 1970. 
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Card Type II--Section Height Card 


Implied 
Card decimal Special 
columns designation Identification notes 
1-3 XXX Plot number 
4 Xx Tree number 
5-7 XXX Species See ADP code table 2, 
appendix I 
8-11 XXX, X D.b.h. (0.b.) Inches 
12-15 XXXX Breast height age Years 
16-17 XX Card number Card position in 
tree deck 
18-19 XX Number of sections (Including tip) 
limit: 30 or less 
20 Xx Ring count measurement 1 = 10-yr interval 
interval 5 = 5-yr interval 
3 = 1-yr interval 
21-24 XXX, X Stump height (ft) Section 1 
25-28 XXX, X Height above the stump Section 2 
of breast height section 
29-32 XXX. X Height above the stump Section 3 
of section 3 
33-36 XXX. X Height above the stump Section 4 
of section 4 
77-80 XXX. X Height above the stump Section 15 


of section 15 


Up to 15 section heights can be punched on one computer card, and the program can 
process up to 30 section heights for any one tree. The additional section heights are 
punched on a continuation card (Same format as above) with section 16 height data entered 
in columns 21-24. It should be noted that the tip height measurement is considered one 
section. Columns 1 to 20 of the continuation card will be a duplicate of the previous card 
except for columns 16 and 17 which will contain the card's position number in the tree deck. 
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Card Type III--Section Diameter (Outside Bark) Card 


For every section height, a corresponding diameter outside bark is needed (tip height 
diameter equals zero). 


Implied 

Card decimal Special 
columns designation Identification notes 

1-15 (Same as section height card) 
16-17 XX Card number Card position in 

tree deck 
18-20 (Same as section height card) 
21-24 XX. XX D.o.b. of stump section Section 1 
25-28 XX. XX D.o.b. of breast height Section 2 
section 

29-32 XX. XX D.o.b. of section 3 Section 3 
33-36 XX. XX D.o.b. of section 4 Section 4 
77-80 XX, XX D.o.b. of section 15 Section 15 


If the number of section heights (including tip height) exceeds 15, a continuation card 
will be necessary with the diameter outside bark of the 16th section entered in columns 
21-24. The first 20 columns of the continuation card will be identical to the previous 
card except columns 16 and 17 will contain the card's position number in the tree deck. 
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Card Type IV--Single Bark Thickness Card 


For every section height a corresponding single bark-thickness is needed (tip height 
bark thickness equals zero). 


Implied 
Card decimal Special 
columns designation Identification notes 
1-15 (Same as section height card) 
16-17 XX Card number Card position in 
tree deck 
18-20 (Same as section height card) 
21-24 XX, XX Single bark thickness of Section 1 
stump section 
25-28 XX. XX Single bark thickness of Section 2 
breast height section 
29-32 XX. XX Single bark thickness of Section 3 
section 3 
39596 XOX OK Single bark thickness of Section 4 
section 4 
77-80 XX. XX Single bark thickness of Section 15 
section 15 


If the number of section heights (including tip height) exceeds 15, a continuation card 
will be necessary with the single bark thickness of the 16th section entered in columns 
21-24. The first 20 columns of the continuation card will be identical to the previous 
card except columns 16 and 17 will contain the card's position number in the tree deck. 
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Card Type V--Section Ring Count 


For every section height a corresponding section ring count is needed (tip height ring 
count equals zero). 


Implied 

Card decimal Special 
columns designation Identification notes 

1-15 (Same as section height card) 
16-17 XX Card number Card position in 

tree deck 
18-20 (Same as section height card) 
21-24 XXXXK Ring count of stump section Section 1 
25-28 XXXX Ring count of breast Section 2 
height section 

29-32 XXXXK Ring count of section 3 Section 3 
33-36 XXXXK Ring count of section 4 Section 4 
77-80 XXXXK Ring count of section 15 Section 15 


If the number of section heights (including tip height) exceeds 15, a continuation card 
will be necessary with the ring count of the 16th section entered in columns 21-24. The 
first 20 columns of the continuation card will be identical to the previous card except 
columns 16 and 17 will contain the card's position number in the tree deck. 
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Card Type ViI--Sequential Radial Increment Measurements Card 


The remainder of the tree data deck consists of card type VI. 


card is: 
Card 
columns 


1-15 
16-17 


18-20 
21-24 


25-28 
29-32 


33-36 


77-80 


Implied 
decimal 


designation 


Identification 


(Same as section height card) 


XX 


Card number 


(Same as section height card) 


XX. XX 


XX. XX 


XX.XX 


XX. XX 


XX. XX 


Radial increment measure- 
ments for stump section 
Radial increment measure- 
ments for stump section 
Radial increment measure- 

ments for stump section 3 
Radial increment measure- 
ments for stump section 4 


Radial increment measure- 
ments for stump 
section 15 


The format for this 


Special 
notes 


Card position in 
tree deck 
Section 1 
Section 2 
Section 3 
Section 4 


Section 15 


If the number of section heights (including tip height) exceeds 15, a continuation card 
will be necessary with the radial increment measurement of the 16th section entered in 


columns 21-24. 


The first 20 columns of the continuation card will be identical to the 


previous card, except columns 16 and 17 will contain the card's position number in the 


tree deck. 
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Coding Format Layout Determination 


Before card type VI can be punched, the radial increment data must be recorded on 
stem analysis data forms in a manner similar to that shown in figure 12 of appendix I, 
and as follows: 


il, 


i) 
e 


Calculate the total number of lines needed to record the stump radial measurements. 

Number of lines — Ring count of stump section (years) 

Ring count interval (years) 

If the answer is not a whole number, round the answer up to the next whole number. 
Figure 12 example: 437/10 = 43.7. Next whole number; in this instance, 44 lines. 
Any decimal except .0 requires such rounding. 
Calculate the number of lines needed to record the radial measurements for the 16th 
section (if there is a 16th section). If the 16th section is the tip, then one line is 
required to record the needed zero. Figure 12 example: 93/10 = 9.3. Round up 
to 10 lines. 
Calculate the total number of lines that will be needed for the tree. Total number 
of lines = number of lines used for heights, diameter, age, and bark + number of 
lines needed for stump section + number of lines needed for 16th section. Figure 12 
example: 8 + 44 + 10 = 62 lines. 
In columns 16 and 17 of the stem analysis data coding forms, number the lines 
consecutively up to the total number of lines calculated in item 3. Use additional 
forms if needed. 


- Calculate the line number where the first radial measurement of section 16 will 


be recorded. 
Line number = Total number of lines - [ (number of lines needed for 
section 16) - 1] * 2 (asterisks used in computer coding format layout 
determination formulas mean "multiplied by." ) 
Figure 12 example: 62 - (10-1) * 2 = 44 
Starting on the line calculated in item 5 and skipping every other line, record the 
radial increment measurements for the 16th section in columns 21-24 of the data 
forms. Last entry for section 16 should fall on line 62. 
For all sections above the 16th section, calculate the line number where each 
section's first radial measurement will be recorded on the data forms. This 
calculation is as follows: 
Line number for section 'i'' = Total number of lines - [(number of lines 
needed for section "i'') -1] * 2 
Where: 
The number of lines needed for section "i"' is calculated by the formula 
in item 1. 
Figure 12 example for section 17: Line number = 62 - [(7-1) * 2] = 50 


. Starting on the line number calculated for section ''i'' in item 7, record section 


"i's" radial increment measurement data on every other line in the appropriate 
columns of the data forms. 
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9. Radial increment measurements for the stump section will begin on the line 
immediately following the last section ring count line (card type V) and will be 
recorded on consecutive lines until the first measurement of section 16 is 
encountered in columns 21-24 of the data form. When this occurs, subsequent 
radial measurements will be recorded on every other line until completion of 
the stump measurements. 

10. For any "k" section 2 through 15, the beginning line number can be calculated by 
by the following formula: 

Line number = Total number of lines - (Lines required for 16th section) - (Lines 
required to record section ''k'') + 1: Where the number of lines needed for 
section "k' is calculated by the formula in item 1. 
Figure 12 example -- Section 2: Line number = 62 - 10 -44+1= 9 

-- Section 3: Line number = 62 - 10 - 42 +1=11 

11. Record sections 2 through 15 beginning on the proper line number (as calculated 
in item 10) for the section being measured. Remember that all lines that have 
section 16 recorded on them will be skipped when recording sections 2 through 15. 
Because the above instructions are involved, we recommend that figure 12 of 
appendix I be studied for a full understanding of the proper format for the program 
input cards. 


Program Output 


Program output consists of punched card output as well as printed output. An example 
of printed output is shown in figures 16 and 17 of this appendix. 


Two card decks are produced for each tree. One deck contains the rescaled sequential 
radial increment data and the other contains the age-height-site index information. The 


rescaled sequential radial increment deck has a format identical to the input card format. 


The format for the age-height-site index output card is: 


Implied 
Card decimal Special 
columns designation Identification notes 
1-2 XX Study 
3-6 XXXX Project 
7-10 XXXX Species code See table 2, appendix I 
11-13 XXX Plot number 
14-16 XXX Tree number 
17-18 XX Composite tree (normally 1) 1 = single tree, 
2 or greater = 
composite tree used 
for averaging 
19-22 XXXX Stump age 
23-26 XXXX Total breast height age 
27-30 XXXX Site index at b.h., age 100 Blank if tree is less 


than 100 years at b.h. 
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Card 
columns 


31-34 


35-38 


39-46 


47-49 


50-80 


Implied 
decimal 


designation 


XXXX 


XXXX 


XXXXXXXX 


Special 


Identification notes 


Height above breast height Interpolated values 
for each decadal 
ring count interval 
multiple 

Cumulative breast height Decadal multiples of 

age by successive decades ring count measure- 
ment interval 

Blank 

Number of cards in this 

tree's site index deck 
Blank 
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PROGRAM NOBLE (INPUT ,0UTPUT,PUNCH, TAPE8,TAPES=INPUT ,TAPE65=OUTPUT, 
1TAPE7=PUNCH) 

DIMENSION HT(31),DIAC31) BI (31) yIAGE (31), RADII (805318) yASAS(31 Dy 
PAGABH(31) ,HTABH(31) sHTTYI (80), I1P (31) yCRAD( 80) yRESCL (31,80) ,IRL(3SI, 
280) »KHT (80) 5 JH( 80) 5 JO( 80) » UB( BD) 

SgKAXKIC3) pXAK2C3) gp VYAXI(03) ,YAX2(03) BUF (2000), 
&CHT (31) 

DATA YAX! /3GH HEIGHT OF SECTION (FEET) 

DATA YAX2 /30H TOTAL HEIGHT OF TREE (FEET) 

DATA XAX! 30H RADII MEASUREMENTS CINCHES) 
DATA XAX2 /30H AGE AT STUMP (YEARS) 

ISKIP=0 

IFCISKIPEQI) GO TO 79% 

CALL PLOTS(BUF,2000, 38) 

CALL SYMBOL (BoepDe go Dg 390 oDge ld 
CALL FLOT (2oe9le 92-3) 

CALL ZSUPRS 

CALL FZSUPR 
ISKIP=I 
CONTINUE 


~ BABA 


READ *°*00° CARD 


READ (Sy IS72VIHP,THT gy IHSy ITHGBH, JXy IHC yg IT HSEC,y THAI, YCUT, 52L~9Z1ly 22523 
19Z%9 25 9269279285 ELEV g ASP yHSLOP gNLAT  NLONGyPNFSI pNERHT gp ISOCRDy 
Z2ISFNT » ISSPT y»IHASRDg IHAPNT,IHAPT 

FORMAT (IS Tt pl 3Sy2l Gg 2l 2g! gF 2c Dg F2eDy &lASy Ald gF5eDyF2eDgF SeDy2IT by 
1A3ys1T4_5Xy,6T1) 

IF (EOF 95)9000,88 

CONTINUE 

YYCUT=1900.*YCUT 


ZERO ARRAYS 


DO 134 NAO=!1,80 

DO 136 NAP=1,31 

HT (NAP) =GIACNAP) =BT (NAP) =AGAS (NAP) =AGABH (NAP) =HTABH(NAP)=0. 
IAGE(NAPD=IFP (NAP) =0 

RADITCNAGgNAPI=RESCL (NAP NAO) =020 

IRL (NAP, NAO) =0 

HTTYI (NAO) =CRAD(NAO) =0.0 

KHT (NAO) =JH (NAO) =JOCNAO) =JBC( NAO) =9 

CONTINUE 

CONTINUE 


READ HEIGHTS OF SECTIONS 


READ(S,1 OODDIIPLIT,ITREE,ISFy,DBHyTAgISECy IAL, (HT (I) yp I= 915) 
FORMAT CIS, Tig lS, Fu old slug 2Xgl2n Tig I SFuel) 

IF CISEC. GT 215) 253 

READ(5,1001) C(HT(I) pI =16,390) 

FORMAT (20X11 5F 4.1) 


oan 


(oil 7 a) 


OanD 


oll a) 


1004 


139 
138 
12 


READ DIAMETERS OUTSIDE AARK OF SECTIONS 


READ(5,1002) (DIACJ) ,J=1,15) 
FORMAT (20Xy915F 42) 
IF(ISEC GT 15) 495 
READ (5,1002) (DIA(J) 9 J=10,30) 


READ BARK THICKNESS (SINGLE BARK THICKNESS) 


READ(S,1002) (BICK) »K=1 515) 
IF (ISEC.GT.15) 5597 
READ(5,1002) (BT (CK) »K=1 6, 30) 


READ SECTION RING COUNT 


READ(S,1003) (IAGE(L) pL =1 915) 
FORMAT (20Xy15I4) 

IF (ISEC. GT. 15) 8,9 
READ(5,1003) CIAGE(L) yi =16, 30) 
IF (IAI eEQ.1)XINC2=10. 

IF (LAI 6EQ.5) XINCR=5S.6 

IF CIAI EQ. 3) XINCR=Ie 
XAGE=IAGE(1) 


READ MEASURED RADIT OF SECTIONS 


JAC=KAGE/XINCR 

EXST=XAGE=(JAGFKINCR) 

IF CEXST. GT. 0) JAC=JACHI 
JAB=IAGE C16) /KINCR 
EXTRA=IAGE (16) -(JAB*XINGR) 

IF CEXTRA GT. 9) JAB=JABEI 

NNON=JAC=JAB 

DO 10 M=1,80 

READ (55,1004) CHEGKy (RADII (MyN) gN=I 515) 
FORMAT CISXyF 2005 3Xy15F 4.2) 

IF (CHECK.E£Q.99.) GO TD 100 
IFCISEC.GT015) 31,1G 

IF (MeGTeNNONI7!1, 10 

READ(5 51004) CHECK, (RADII (49N),N=16,30) 
IF (CHECK.EQ.99-) SO TO 100 

CONTINUE 


CALCULATE NUMBER OF CARDS 


CARDS=XAGE/XINCR 
ICARD=CARDS 

XPAGc=IAGE(2) 
CROS=XPAGE/XINCR 
IGRDS=CRDS 

IF (CRODS-ICROS) 138,138,139 
ICFCS=1ICPOS+! 

IF (CAxXDS-ICARDI 13,1312 
ICAxD=iCARD+1 


45 


Haan 


OO 


3a 


14 


784 
782 


783 


20 


19 


17 


1402 


1603 


46 


CALCULATE AGE JF TREE AT A GIVEN HEIGHT ABOVE STUMP 
AND S8KEAST HEIGHT 


DO Ie IJ=iygISEC 

TI J=IJ+! 

AGAS(IJ) =LAGEC I) -TAGE(I SJ) 

AGABH(LIJI=LAGE (2) -IAGE (II) 

CALCULATE HEIGHT ABOVE S3SREAST HEIGHT 

HTABH(IJ)=HT (IIS) -HT (2) 

CONTINUE 

CALCULATE HEIGHTS ABOVE BREAST HEIGHT AT SELECTED YEAR INTERVALS, 

NOSEC=] 

TVY=XINCR#.95 

NX=IAGE(2)/1Y+#1 

AXL=NX*I¥= ITAGE(2) 

TF CAXL -EQ-LTY) NX=NX=1 
DO 15 IK=tyNX 

IF(IK.EQeNX) GO TO 16 
COAGE=XINCR*IK 

IF (COAGE oGT -AGABSH(NOSEC)) 782,783 
NOSEC=NOSEC#) 

GO TO 784 
HTINT=HTASH (NOSES) -HTABH(NOSEC =I) 

AGINT=AGABH(NOS=C) -AGABH(NOSEC=1) 
DIF=COAGE-AGABHI(NOSEC=1) 

HT ADD=(DIFSAGINT) *HTINT 
HTTYICIK)=HTAPRH(NOSEC-1) +HTADD 

GO TO Is 
HT TYICIK)=HTABH(ISEC=1) 

KHT (IK) =HTTYICIK) +25 
CONTINUE 
CALCULATE RESCALE OF RADII FOR INDIVIOUAL SECTIONS 
0O !7 ILL=!,ISES 

IFCILLeEQeISEC) GO TO 17 
BO 19 IKO=!1,ICARD 

IC=ICARD 

IP (ILL) =ICARD 
CRADCILLI=(OIACIL LI -2.*8TCILLII/2. 

RESCLCUILL »pIXDI=H(CRAN(ILLISPADITICIC,ILL)) *XADII(IXD,IL) 

IRL (ILL » IXO) =RESCL AIL gs IXDOI*I00.+.5 
CONTINUE 

IF(IHAPNT.GT.0950 TO 1608 

WRITE(6,1102) 

FORMAT (1H! 954Xy3HOUMS24XKX,35HC STUMP 8eHe SITE HEIGHT BeHe / 
L734 ST PRO SPES CPEOT Oi Rec. Ta eAGE AGE 405 ~4ed AGE 
(2 CARDS 4 ) 

CONTINUE 


DO 23 KIX=1,I1C 
MAGE=KIX*X INCRE 
IF (MAGE. GT.ITAGE (2) ) MAGE=IAGE C2) 


aOaqa 


1100 
1602 


PUNCH OUT SITE INDEX CARD 


NS=0 

SKKHER=HTTYICI0) 

IF CIAGEC2]) -LT2100)SKKHER=0. 
KKHER=KHT (40) 

IF CIAGE(2) -L T2100) KKHER=0 

IF (IHACRO.GT.0)GO TO i602 

WRITEC7 yA 1OOYNSyISPyITPLOT,ITREEsITAGE (ID yIAGE(2) »KKHERy HT (KIX) 
IMAGE,ICRXDS 

FORMAT (2H 1y2LT4y2T3y2H 1514 98X_13) 
CONTINUE 

IFCIHAPNTeGT.09GO TO 1603 


$1015 WRITECG,1ITONINS sISP,IPLOT,ITREE,TAGE (I) y TAGE (2) ySKKHER, HT TYICKI 


oan 


1103 
(603 


23 
22 


1103 


1X) yMAGE,ICROS 

FORNAT (2H 19415, 3H Fy l5 yl 6y2FB.2 9169 8XyI5) 
CONTINUE 

IF (MAGE*EQeIAGE(2)) GO TO 22 

CONTINUE 

CONTINUE 


WRITE AND PUNGH OUT INFORMATION ON RESCAL= 


HERHT=NERHT* ol 

HLONG=NLONG*.! 

HLAT=NLAT*.1 

DHOBH=IHDBH*.4J 

ZLONG=NLONG*.003 

ZLAT=NLAT*.001 

JLON=ZLONG 

JLAT=ZLAT 

ZLONG=JLON* 100.0 

ZLAT=JLAT*100.0 

ZLATM=HLAT=ZLAT 

ZL ONGM=HLONG=ZLONG 

HLAT=JLAT 

HL ONG=JLON 

IFCISFNT2GT.9)G0 TC 160% 

WRITE (06,1573) DHD 3H, SXy CRLL gZ 19 22g Z3yleygPNF Sly 25 yoy Z7gZByHERHT, ELE 
§1VyASPyHSLOP yHLAT »ZLATMyHLONGyZLONGM 


1573 FORMAT CIHI1,42Xy2BHTRES AND PLOT INFORMATION S//5XKy4HD3BH=_9F5ely 


SRELIMINARY NOBLE FIR SIVE INDEX=9A3,10X,I13HCL CAR ~ENGTH=,A3_1He Al 
4yg5KygISHOROWN LENGTH=9A3q1Hoe Al pSXgIBHTOTAL TREE HEIGHT=5"-6.17/ 

SOX gIOHELEVATION=9 5.09 5X g7HASPECTHgF 300g SK gGHSLOPCH=5F Ge Dy IH% 9 SX IS 
OHLATITUDE NORTH=9-4%e094H DEG y= 60e1,23H MIN LONGITUDE WEST=59F 4.05 
74H DEGyF Gel y4H MIN/S) 

WRITE CO,1I 03) ISP,IFPLOT,LTREE (HTC I) g T=1, 21), (DIACI) gy J=1921),08T(K) 
}yK=1,21) 5 (TAGE CKROA) gy KKB=1, 21) 

FORMAT ( D4XsIHSFECIES =5IT4/747XgSHPLOT NOogIT¥aeXyIATXEE NOoyl3 
1/7 7TXy8HSTZHT t+yOXy2HHDISTANCE FROM STUMP,32X,20HDISTAVNCE FROM 
ZSTUMP/7X yFOel sy IXy20F 6015 1X 
3 IX, NS2CIHHISIX,EHDIA OR gy Fe lg I Xg20F Oe IWIXgoHBARK gFOe2gIXy 
420F60271X,1S2(THt)/4H AGE, 3Xy,l6yEX,20I6) 


47 


1604 


4G] 


1105 


42 


43 


4k 


45 


46 


47 


48 


49 


1104 
1605 


48 


GONTINUE 

JOSH=CBHF1I00+.5 

DO 4&t IV=!,ISEC 

JHCIVYI=HHT CIV) *10.¢+.5 

JDCIV) =BIACIV) *1008.+.5 

JBCIVI=BT(CIVI*FIG0.+.5 

KNC= I 

KDBH=D8H* 10.4.5 

IFC(ISCRO.GT.0IGD TO 49 

WRITEC 7, IS72)LHP,ITHT, IHS, LHOBHs JX, IHG, IHSECgIHAL,YCUT »pCRL 9Zl yZ2qZ 
139 Z2&9Z59ZOyZ7_yZB_pELEWy ASP yg HSLOF yNLAT » NLONGyPNFSI 9 NERHT 

WRITEC7, SV OSVIPLOT pITREEgISPyKDSHyIAg KNGgISEC,IALs (JHC KV) y KVH1515) 
FORMAT (13, Ti 913,214, 21251!) 51514) 

LFCISEC.GT.35)G9 TO 42 

GO TO 4&3 

KNC=KNC#I 
WRITEC7gLIOSIIP.LOT pI TREE gI SF a KDBHy ITA, KNUgISECsTIALy (JHC KV) y KV=H16, 30 
1) 

KNC=KNC#] 

WRITEC 7, IV OSI IP_LOT pS ITREEy ISP y KOSH ,ITAgKNU,ISEC yi AL, (JOC KW) y KWH=1, 15) 
IF (ISEC.GT~15) GO TO 4&s& 

GO TO 45 

KNC=KNC#I 

WRITEC7, 1105) IP OT pI TREE I SF y»KDBH,IA,KNC,ISEC sIALs (UDI KW) yKW=16, 30 
1) 

KNC=KNC#1 

WRITE(7 JI DSIIPLOT pLTREE oI SP y KDBH TA, KN yg TSEC gTALy (CIB KX) 5 KX=1515) 
IFCISEC.GT.15) GO TO 46 

GO TO 47 

KNC=KNC# | 

WRITE(7, JIOSVIPLOT pI TREe »yITSPy»KDBHyIAgKNCO gL SEC IAL y (UB KX) gKX=165 30 
1) 

KNC=KNC #4 

WRITEC7, 1195) IPLOT ,ITREEs ISP yKDBHyTAgKNGsIScOsgIALs (IAGE(KS) pKS=ly 
415) 

IFCISEC.GT~15) GO TO 48 

GO TO 4&9 

KNC=KNC#I 

WRITEC7 INOS) IPLOTsI TREE gI SP gKDBHyITAgKNGgISEC TAI, (IT AGE(KS) yKS=16,y 
139) 

TADI=IAGECI) SXINCR 

LOV=IAGE (1) -(XINCPFIADI) 

IF (LOV.EQ.0)IADI=IADI-$ 

ITAD2=I ADI 

DO 2& IPOJ=!1,ICARD 

IHIST=VYCUT=(TADI*XINCR) 

IADI=IADI-I 

IFCISFENT.GT20)G9 TO 1005 

WRITE (6, 1104) IHIST y(RESCLCIByIPOJ) ,IB=1, 21) 

FORMAT (IS) 2X pF 60291X g2DF 60 2) 

CONTINUE 

KNC=KNC#t] 

IF (KNC -GT.98) KNO=986 

IF(ISCRO.GT.VIGD TO 24 

WRITEC7T, I IOSVIPLOT,ITREc sISPyKDBH yg IAgKNCgISECsIAL, (IRL (IG,1 P04) IG 
1=1,15) 


51 


52 


24 


53 


1107 


HEIGL S ECs Glico) GOr TO. 5H 

GO TO 24 

IFCIRL (t6,TPOU)-GT.0) GU TO 52 
GO TO 24 

KNC=KNC#1 

IF (KNC oGT. 98d KNC=58 


WRITEC7y 11 O5) IPLOT,ITREEgISPyXOBHyIAgKNC,ISECyIALT, (IRL (15,1P05) ,IG 


=16,30) 

CONTINUE 

IF(ISES.GT.21) GO TO 53 
GO TO 8787 

CONTINUE 

IF CISPENT.GTe0)G9 TC 3787 


WRITE CG, IIO7IISP,IFLOTgITREEg (HT (1) yg LTH22930) 5 (DLAC) y J=22y5 30) g (BT 


1K) sK=22, 30) 5 (TAGE (KRB) »KRB=22, 30) 


FORMAT CTH! s53X,9HSPECIES = 9I4/47Xg9HPLOT NOcgl4a%X59HTREE NOoygl3 


194X s39HCONTINUED//I3X,21HDISTANCE FROM STUMP,32X,2I1HDISTANCE 
2M STUMP/ISKy9IFG0eI1SIX, SO7CIH*H)/SIXZE6HDIA OB56XyS9FH0I SIX, H04BARK 
Se9FGe2FIX,6701HtIS4H AGE9X, 916) 


DO 54% IU=!,ICAR) 

ITHIST=YYCUT-(TAD2* xXINCR) 

ITAD2=IAD2-1 

IF (RESCL (22,1U)-£Q.0.9) GO TO 54 
WRITE(6,/108) THIST gy (RESCL(JUgIU) g JU=22, 30) 
FORMAT (1558X59F 5.2) 

CONTINUE 

CONTINUE 

IFCISSPT.GT.0)GO TO 1609 

YINK=(0. 

IF (IAL.EQ.3) VINK=5.0 

XINK=2.0 

IF (IAI~EQ.3) XINK=1.0 

CALL AKIS(0.0,0.09VAXKI1 93092705590 20,0.0,VYINK,I10.5-1) 
CALL AXIS(0.0,0e0yXAXI y-350520.60,0060y Oe DgXINKy2D05-1) 
SPI=ISP 

Tine ek ihe = 

PLOTI=IPLOT 

SSTAGE=I1 AGE (1) 

CALL SYMBOL (5.0527 00902551 2HSTEM PROFILE 90.0512) 
CALL SYMBOL (500g 260590 2594HPLOT 9 0005 4%) 

CALL NUMBER (6025926059 e259PLOTI_g0.0,-14) 

CALL SYMBOL (5.05 20-0592 2554HTREE 90.05%) 

CALL NUMBER (6025926005 025, TREEI yO. 0y=-1) 

CALL SYMBOL (52.039 254592 2597HSPECIES 900097) 

CALL NUMBER (7.255925 055 e259SFI9 0005-1) 

CALL SYMBOL (5405254092255 3HAGE 5 020,53) 

CALL NUMBER (6625925005 025ySSTAGE,000g-1) 

CALL SYMBOL (5.0924.55425y6HBH AGE,0.090) 

CALL NUMBER (70259 2405590 25 9XPAGE9 0009-4) 


49 


FRO 
9 OX 


aaonn 


OO 


50 


447 


642 


644 


449 


448 


452 


CALCULATE HEIGHT BETWEEN SECTIONS FOR STEM PROFILES 
AND FLOT STEM ?ROFILE 


DO 445 IHER=1,80 

DO 444 IREX=1,31 

CHT (CIREX) =HT (1) tHT CIREX) 

IF CRESCL (ly IHERD -EQ-0.0) GO TO 455 
IF(IRecX.tQ.1) CHT CII=HT (I) 
AXISH=CHT(IREXI/VINK 

IF CRESCL (TREX, THER) LE 0000) 447 5448 
DIFFA=ITAGECTREXK=1) -TAGECIREX) 
OIFFH=CHTCIREX) -CHTCIREX=1) 
ISA=IAGE (I) /XINGR 
DA=ITAGe(1)-CISA*XINCR) 

IF (DAcEQ e000) 641 9542 
SA=DA+4+(THER-1)*XINCR 

GO TO 449 

SA=(IHER) *XINCR 

ANOM=SA-AGAS(IREX=1) 

AXISH= (C(CANOM/DIFFA) FOITFFH) tCHT (TREX-1))/7YINK 
AXISR=0.0 

CALL SYMBOL CAXTSRy AXISHy 0 059350005-2) 
GO TO 445 

IF CAXISHeGTe27e5) GO TO 455 
AXISR=RESCL (CIREX gy THER) /XINK 

IF (CAXISReGT.2020) GO TO &55 

IF CIREXeEQe1) 4595452 

GO TO 444 

CALL SYMBOL (CAXISR ey AXISH, 005935000 9-2) 
CONTINUE 

CONTINUE 


PLOT OUT HEIGHT OVER AGE GRAPHS 


CONTINUE 

CALL ADVANZ(44) 

CALL PLOT(00050.0519) 
CONTINUE 
IFCIHAPT.GT.9)G0 TO 1 
YINK=10.6 

IF (IAT .EQe 3) VINK=5.0 
XINK=25.20 

IF (IAI .EC.3) XINK=5.0 


CALL AXIS(DeDyBeDsVAX2, 30g 270599000, 00DgVYINKgIOD.0 5-1) 


CALL AXTS(0 eI 9 De Dg XAX2g- 3D g2100y DeDg DoD yXINKyI0009-1) 
CALL SYMBOe (50 Dy 2605 ye 25, 4HPLOTy D0 0g 4) 

CALL NUMBER (6025920057025 gFLOTI 9 V0D5-1) 

CALL SYMBOL (5009250092259 4HTREE9 Do Dy &) 

CALL NUMBER (6025526009 e25gTREELTs 000g- I) 

CALL SYMBOL (5605 25059025 y7HSPECIES yO 0097) 

GALL NUMBEX (70255925055 eo 25ySPIy 00 Dy =I) 

CALL SYMPOL (5209250090 259 3HAGE » De Dy 3) 

CALL NUMBER (6025925009 e2SySSTAGE 000g - 1) 


453 
456 


9000 


986 


CALL SYMBOL (5209240556259 6HBH AGEy060y5) 


CALL NUMBER (7 e252 ed ye 2dgXPAGE 0205-1) 
DO 453 IREG=1!,31 

XP=AGAS(IP2EG)/XINK 

YF=CHT (ILREGI/SYINK 

IF (AGASCIREG) -GT2£525) GO TO 456 

IF (CIAGE(IREG) eLE©OoAND.HTCIREG) LED. D0) 
IF (XP eGT.212.0) GO TO 456 

IF (YP.GT.28.-0) GO TO 4506 

CALL SYMBOL (XP YP, 20535906 0y-1) 
CONTINUE 

CALL ADVANZ(44) 

CALL FLOT(06050.0,19) 

GO TO ! 

CONTINUE 

CALL PLOT (2220590.05-3) 

CALL PLOT(0,0,999) 

STOP 

END 


GU 


Tr 
' 


0 


456 


ol 
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The mission of the PACIFIC NORTHWEST FOREST 
AND RANGE EXPERIMENT STATION is to provide the 
knowledge, technology, and alternatives for present and 
future protection, management, and use of forest, range, and 
related environments. 


Within this overall mission, the Station conducts and 
stimulates research to facilitate and to accelerate progress 
toward the following goals: 


1. Providing safe and efficient technology for inventory, 
protection, and use of resources. 


2. Developing and evaluating alternative methods and 
levels of resource management. 


3. Achieving optimum sustained resource productivity 
consistent with maintaining a high quality forest 
environment. 


The area of research encompasses Oregon, Washington, 
Alaska, and, in some cases, California, Hawaii, the Western 
States, and the Nation. Results of the research are made 
available promptly. Project headquarters are at: 


Fairbanks, Alaska Portland, Oregon 
Juneau, Alaska Olympia, Washington 
Bend, Oregon Seattle, Washington 
Corvallis, Oregon Wenatchee, Washington 


La Grande, Oregon 


Mailing address: Pacific Northwest Forest and Range 
Experiment Station 
P.O. Box 3141 
Portland, Oregon 97208 


GPO 994-595 
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